Mon. Not. R. Astron. Soc. OOP. [THTSl (2011') Printed 31 October 2011 (MN WI^ style file v2.2) 

The Connection Between Thermal and Non-Thermal 
Emission in Gamma-ray Bursts: General considerations 
. and GRB090902B as a Case Study 
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ABSTRACT 

Photospheric (thermal) emission is inherent to the gamma-ray burst (GRB) "fireball" 
model. We show here, that inclusion of this component in the analysis of the GRB 
prompt emission phase naturally explains some of the prompt GRB spectra seen by 
the Fermi satellite over its entire energy band. The sub-MeV peak is explained as 
multi-color black body emission, and the high energy tail, extending up to the GeV 
band, results from roughly similar contributions of synchrotron emission, synchrotron 
self Compton (SSC) and Comptonization of the thermal photons by energetic elec- 
trons originating after dissipation of the kinetic energy above the photosphere. We 
show how this analysis method results in a complete, self consistent picture of the 
physical conditions at both emission sites of the thermal and non-thermal radiation. 
We study the connection between the thermal and non-thermal parts of the spectrum, 
and show how the values of the free model parameters are deduced from the data. We 
demonstrate our analysis method on GRB090902B: We deduce a Lorentz factor in the 
range 920 ^ 77 ^ 1070, photospheric radius r^/j ~ 7.2 — 8.4 x 10^^ cm and dissipation 
radius r-^ ^ 3.5 — 4.1 x lO^'^ cm. By comparison to afterglow data, we deduce that 
a large fraction, e^ « 85% — 95% of the kinetic energy is dissipated, and that large 
fraction, '^ equipartition of this energy is carried by the electrons and the magnetic 
field. This high value of e^ questions the "internal shock" scenario as the main energy 
dissipation mechanism for this GRB. 

Key words: gamma raysibursts — gamma rays:theory — plasmas — radiation 

mechanisms:thermal — radiative transfer — scattering 



1 INTRODUCTION 

Although extensively studied for nearly two decades, the ori- 
gin of the prompt emission of gamma-ray bursts (GRBs) is 
still puzzling. Up until these days, GRB prompt emission 
spectra are often modelled as a smoothly broke n power-law, 
which is referred t o as the "Ba nd" function (JBand et al. 
19931: iPreece e^ al. 1998a, 2000,: iKaneko et al. ■ 2006, 2008; 



Abdo et aL\\2009\i ). While the "Band" function often pro- 



vides very good fits to the observed spectra over a limited en- 
ergy range, it suffers two crucial drawbacks. In several bursts 
seen by the Fermi satellite, a high energy tail extending up to 



tens o f GeV was seen (e. g., GRB0905 10, Ackerman et al. I 
(|2010t ): or GRB090902B. lAbdo et oTI (|2009al ')'). The first 
drawback is that this tail requires more than the "Band" 
function on its own to have an acceptable fit. However, the 
most severe drawback is that the "Band" function, being 
empirical in nature, does not provide any information about 
the emission mechanism that produces the prompt radia- 
tion. 

A common interpretation of the observed GRB 
spectrum is that it results from synchrotron emis- 
sion, accompanied perhaps by s ynchrotron-self Compton 
(SSC) emission at high energies (jRees fc Meszaros 1 11994 : 



A. Pe'er et. al. 



Sari fc PirarTlllQQTl: Ipilla fc Loeb 1 1 19981: iGuetta fc Granot I 



2003al : iPe'er fc Waxman I |2004| : (Gupta fc Zhang I l200'i1 ). 
The emission follows the dissipation of a kinetic en- 
ergy. The preval ent dissipation models i nvolve either 
intern al shocks (JRees fc Meszaros I 11994 : ISari fc Piran I 
[l997b|), magnetic energy dis s ipatio n in Poynting dom- 
inated outflows (iThompson I Il994l: ISpruit et al.j I2001I : 



( Thonips 



iGiannios fc Spruit "2005':'Z hang fc Yan Il2011^ ■ or coUisional 
heating (|Beloborodov 2010 '). These models have two main 
advantages. First, they can explain the complex lightcurves 
often seen (albeit with very little predictive power). Second, 
they naturally account for the non-thermal spectrum. 

In spite of these advantages, in recent years evidence 
is accumulating for serious difficulties in these models. A 
well known deficiency of the internal shock scenario is 
the low efficiency of energy conversion (Mochkovich et al. 
1995; Kobayashi et al. 1997; Daignc fc Mochkovitch 1998,; 



Lazzati et 0/71119991 : iGuetta et 0/71120011 : iMaxham fc Zhang I 



20091 ) ■ resulting from the fact that only the energy asso- 
ciated with the differential motion between the expanding 
ejecta shells after the GRB explosion can be dissipated, and 
that this energy is much lower than the energy associated 
with the bulk motion. This is in contrast to the observations 
which show high ef ficiency in 7-ray emission, of the order 
of te ns of percent (jZhang et al. 1 120071 : iNvsewander et al. I 
I2OO9I ). A second drawback is that optically thin synchrotron 
and SSC emission cannot account for the steepness of the 
low eriergy spectral slopes JCrider et al. 1119971 : IPreece et al. I 
Il998bl . I2OO2I ! iGhirlanda et al. I l2003l ) Talthough part of the 
observed steepening may be accounted for when SSC in 
the Klein Nishina limit is considered; see iDaigne et al. I 
l2009l : lBosniak et al. llJOOgh . Third, the high e nergy spectral 
slope varies significantly fr om burst to burst IjPreece et al. I 
Il998al : iKaneko et al. 1120061 ) , with some bursts showing very 
steep high energy spectral slopes. This is in contradiction 
to the expected spectral slope of synchrotron emission from 
a power law distribution of electrons, which is expected to 
produce a fairly flat spectrum, uF^ x u°. Finally, the flt- 
ting is often made to the time integrated s pectrum. An alysis 
of time resolved spectra done by Crider et al. I ([1993) Eind 
IGhirlanda et al. I (2003), has shown that neither the syn- 
chrotron nor the SSC models can explain the time resolved 
low energy spectral slopes. 

These well known difficulties of the synchrotron 
emission model have motivated the study of alterna- 
tive scenarios. T hese include, among other s, quasi-thermal 
Comptonization IIGhisellini fc Celotti |l999|). Compton dra g 



IZdziarski et al. Ill99ll: IShemi I Il994l : iLazzati et al. IIJOOol ). 
jitter radiation fMed vedev M200d ), Compton scattering of 
synchrotron self absorbed photo ns (jPanaitescu fc Meszaros I 
l200d : IStern fc Poutanen 1 12004| ) and synchrotron ern ission 
from a decaying magnetic field (jPe'er fc Zhang Il200q ). 

An alternative model, which is arguably the most 
natural on e, is a radiative contribu ti on from the pho- 
tosph er e IIEichler fc Levinson I2OO0I: 



200c ; iMeszaros et al. 



200Z: Rees 



c Meszaros 
2010; loka 



12002 

'2005'; 

2010; 



iMeszaros fc Rees 



Daigne fc Mochkovitch . 
Ghirlanda et al. I ^OOtT 



Beloborodov 2010; loka 2010; Mizuta et aQ I2OIOI '). In- 
deed, a photospheric emission is a natural outcome of the 
"fireball" model. At small radii, close to the inner engine, 
the o ptical depth is huge, r ~ 10^^ (see, e.g., IPiran I 
I2OO5I . for a review). As a result, photons cannot escape. 



but are advected with the fiow until they decouple at the 
photosphere. The huge value of the optical depth implies 
that regardless of the initial emitted spectrum, if photons 
are emitted deep enough in the flow, the spectrum emerging 
at the photosphere is black-body (a Wien spectrum may be 
obtained if the number of photons is conserved). In fact, in 
the classical "fireball" model, photons serve as mediators for 
energy conversion (the explosion energy is converted into 
kinetic energy of the relativistically expanding plasma jet 
by scattering between photons and leptons, which in turn 
convert the energy to the baryons via Coulomb collisions). 
Therefore, the appearance of thermal photons as part of 
the observed spectrum is not only expected, but is in fact 
required. 

Apart from being an inherent part of the "fireball" 
model, one of the key advantages of the idea that at least 
part of the observed spectrum originates from the photo- 
sphere is that it provides a relief for the effici e ncy problem of 
the i n ternal shock raodel JPe'er et al. |[2005l : iRvde fc Pe'er I 
I2OO9I : iLazzati et al. I [2009|) . This is because if indeed part 
of the photons that we see have photospheric origin, than 
the total energy seen in photons is higher than the energy 
released by the internal dissipation: only the non-thermal 
part of the spectrum is required to originate from energy 
dissipation above the photosphere. Thus, the dissipated en- 
ergy may be smaller than the remaining kinetic energy, even 
though the total energy seen in the photon component is 
comparable to the kinetic energy. A second great advantage 
of this idea is that the Rayleigh- Jeans part of the thermal 
spectrum (or modification of it) can naturally account for 
low energy spectral slopes much steeper than those allowed 
by the synchrotron theory, and are hence consistent with 
observations. 

The observed GRB spectra are therefore expected to 
be hybrid, i.e., containing both thermal and non-thermal 
parts. These two parts are connected in a non-trivial way. 
As the non-thermal part originates from energy dissipation 
(one or more) occurring above the photosphere, it is natu- 
rally delayed with respect to the thermal photons. A com- 
plete decomposition of the thermal and non-thermal parts 
of the spectrum is in fact nearly impossible. Thermal pho- 
tons serve as seed photons to Compton scattering by the 
energetic electrons, thereby affecting the non-thermal part 
as well. Even if the dissipation occurs high above the pho- 
tosphere, thermal photo ns significant l y con tribute to the 
cooling of the electrons IjPe'er et al. 1 12005| ). therefore af- 
fecting t he spectrum emitte d by these electrons. As was 
shown bv lPe'er et al. I (|200q ). under these conditions a fiat 



energy spectrum (vF^, 



, resulting from multiple Comp- 



ton scattering of the thermal photons, is naturally obtained 
for a large range of parameters. A similar conclusion was 
drawn in the case of a magnetized outfiow under the as- 
sumptio n of slow heating due to continuous reconnection at 
all radii (Gian nios 1120061 ). These predictions were shown to 
be in qualitatively very good agreement with Fermi results 
(|Toma ei al. II2OI0I ). 

The non-trivial connection between the thermal and 
non-thermal components, combined with the fact that both 
components vary in time, make it difficult to identify the 
thermal component. A breakthrough in identifying this com- 
ponent was made bv lRvde 1 1|2004| . |2005| ). who looked at time 
resolved spectra, thereby allowing an identification of the 



Connection between thermal and non-thermal emission in GRBs 



temporal evolution of the temperature. Repeating a simi- 
lar anal}^sis_on_afairljHxirge sample of bursts, it was shown 
by iRvde fc Pe'er I (2009) that the thermal component not 
only could be identified, but that both the temperature and 
flux of this component show repetitive temporal behavior at 
late times: Egg oc i"^, and T°^ oc t"^/^. This behavior was 
found to be in very go od agreement with the theoretical pre- 
dictions (jPe'er 1120081 ) . thereby providing an independent in- 
dication for the correct identification of the thermal emission 
component. It should be pointed out here, that the thermal 
emission is expected to appear as gray-body emission, com- 
posed of multi-color black body spectra. This results from 
the fact that at a given time interval, an observer sees pho- 
tons emitted from different radii and different angles, hence 
undergoing different Doppler shifts. A full analysis shows 
that the resulting low energy spectrum (below the observed 
spectral peak, typica lly seen at sub MeV ) is a power law 
(|Pe'er fc Rvdell201ll ). 

Once the thermal component is identified, it is relatively 
easy to use its properties to deduce the dynamics of the out- 
flow. As opposed to the non-thermal part, whose emission 
radius is uncertain (there may be multiple emission radii), 
the emission radius of the thermal photons is defined to be at 
the photosphere. Thus, by studying the properties (thermal 
flux and temperature) of the thermal component, under the 
assumption of constant outflow velocity, it is possible to de- 
duce the photospheric r adius and the Lorentz factor of the 
flow l|Pe'er et al. II2007I ). We note that this method is in- 
dependent, and is complementary to the opacity argument 
method commonly used to constrain the Lorentz factor of 
the flow at the emission radius of the high-energy photons 



dSvensson lll987l:]Krolik fc Pier Ill99ll : I Woods fc Loebll 19951 : 



iLithwick fc Sari II2OOII ). It has two main advantages to the 
opacity argument method. First, it provides a direct mea- 
surement of the Lorentz factor, rather than a lower limit. 
Second, it is independent of measurements of the variability 
time, which is highly uncertain. 

While in the past a clear identiflcation of a thermal com- 
ponent was difficult, the situation dramatically changed with 
the broad band spectral coverage enabled with the launch of 
Fermi. Out of 14 GRBs detected by the LAT until January 
2010, 3 show clear evidence for a distinctive high energy 
(> MeV and up to the GeV range) spectral component, 
and another 8 show marginal evidence for such a component 
(iGranot 2010 ). Thus, it is natural to deduce that the low 
energy photons (below and at the sub MeV peak of the 
ffux) have a different origin than the higher energy (LAT) 
photons. 

Out of the 3 bursts that show clear evidence for a dis- 
tinct high energy spectral component, GRB090902B may be 
the easiest to analyse, for two reasons. First, the low energy 
part of the spectrum (the spectral peak) is very narrow, 
and both the low energy and high energy spectral slopes are 
very steep. Thus, any attempts of fitting the spectrum us- 
ing the standard synchrotron-SSC model are rejected. On 
the other hand, the spectrum is easily fitted with a (multi- 
color) black body spe ctrum, plus an additional power law 
l|Rvde et al. I I2OIO: Zhang^rfliDlioil). Thus, to date, this 
burst is unique by the fact that a thermal component is so 
clearly pronounced in its spectrum. Second, the high energy 
power law component is spectrally distinctive from the low 
energy one. Using the opacity argument, the detection of 



a 33.4 GeV photon associated with this burst (JAbdo et al. I 
l2009al ) implies that the emission radius of this photon must 
be much greater than the photosphere (see below). These 
two facts make this burst ideal for demonstrating how sep- 
aration of the spectrum into thermal and non-thermal com- 
ponents enables one to deduce the physical conditions of 
the outflow at the emission sites of both the thermal and 
non-thermal photons. Moreover, as we will show below, one 
can use measurements of the non-thermal part, to remove 
some of the uncertainties that exist in measurements of the 
thermal component alone. 

This paper is organized as follows. In section 321 we 
provide a general discussion on the properties of the ther- 
mal - non thermal model. We show how one can combine the 
hydrodynamic information derived by studying the proper- 
ties of the thermal component with the constraints given 
by measuring the properties of the non-thermal part of the 
spectrum, to provide a comprehensive picture of the physi- 
cal parameters at both emission sites. In !j3]we demonstrate 
the use of our method by fitting the prompt emission spec- 
trum of GRB090902B, and deducing its physical properties. 
We summarize and conclude in ^ 



2 GENERAL PROPERTIES OF A THERMAL - 
NON-THERMAL MODEL 

2.1 Temperature, luminosity and constraints on 
the value of te 

We consider a fireball wind of total luminosity L, expand- 
ing from an initial radius ro (which, for the sake of argu- 
ment, can be assumed to be a few times the last stable orbit 
around the central black hole, or the sonic radius; in any 
case, it is a few times the Schwarzschield radius around a 
non-rotating black hole; see further discussion at the end 
of section [5TT1) . The initial black-body temperature at ro is 
To — {L/ATvrQcay'*, where c is the speed of light and a is the 
radiation constant. As the optically thick (adiabatic) wind 
expands, the baryon bulk Lorentz factor increases as F oc r, 
and the c omoving temperature drops as T' oc r~^ (see, e.g., 
iMeszaros 2006, for a comprehensive review), [j As long as 
the wind remains optically thick, the acceleration continues 
until the plasma reaches the saturation radius, Vs = rjro, 
above which F coasts to a value equal to the dimensionless 
entropy, i] = L/ M(? . Here, M is the mass outflow rate. 

The photospheric radius is the radius above which the 
flow becomes optically thin to scattering by the baryon- 
related electrons. Depending on the values of the free model 
parameters {r],L, and ro), this radius can be above o r be- 
low the saturation radius (see Mcsza ros et al. 1120021 ). For 
the parameter values characterizing GRBs (see below), the 
photosph eric radius is above the saturation ra dius, and is 
given by (JAbramowicz et al. |[l99ll : IPe'er II2OO8I ) 



Tph = 



Lctt 



8'Kmpri'' 



= 5.8 X 10"L54%^cm. 



(1) 



Here and below, or is the Thomson cross section, rup is the 
proton mass, and we use the convention Q — lO^Qx in cgs 



^ From here on, quantities measured in the comoving frame are 
primed, while unprimod quantities are in the observer frame. 
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units. The high values of the luminosity and the Lorentz 
factor chosen for the demonstration in equation [l] are for 
ease of comparison with the Fermi results of GRB090902B 
(see ^ below). 

Above the saturation radius, adiabatic energy losses (in 
the absence of dissipation) cause the temperature to drop 
as T = ro(rs/r)~^". The observed temperature of photons 
emitted at the photosphere is therefore 



yob ^ y^ 



(rph\ 



-2/3 



= 3.6x10^(1+^) 



-1 r-5/12 



% 



8/3^1/6 



eV, (2) 



where z is the redshift. Note the very strong dependence 
of the observed temperature on the asymptotic value of the 
Lorentz factor, r): for high t), high values of the temperature 
are expected. 

The observed photospheric thermal luminosity drops 
above the saturation radius as Lrnir) = {L/2){r/rs)~^'^ , 
the greater p art of the energy being in a kinetic form , L^ ~ 
L/2 (.MeszaroTfc Rees ll200(]| : lRees fc Meszaros l2005l fl. The 
non-thermal part of the spectrum results from dissipation of 
the kinetic energy. Part of the dissipated energy goes into 
accelerating electrons, that radiate the non-thermal spec- 
trum. Denoting by ed the fraction of kinetic energy that is 
dissipated and by te the fraction of dissipated energy that is 
converted to energetic electrons, one obtains an upper Ivmit 
on the ratio of non-thermal to thermal luminosity in the 
spectrum: 






r oh 



EdCf 



m 



2/3 



0.33L5f^73■**/^ 



-2/3 

0,8 '^dfi^e,- 



(3) 



where the non-equality results from the fact that the elec- 
trons do not necessarily radiate 100% of their energy. Equa- 
tion [3] implies an interesting result: the higher the Lorentz 
factor of the flow is, the more pronounced its thermal lu- 
minosity is expected to he. This result may be very signif- 
icant given recent Ferroi-LAT data, which show evidence 
for high values of the Lorentz factors in several bursts. We 
further point out that as rj increases, the saturation radius 
rs increases, while the photospheric radius decreases. Thus, 



at high enough value oi rj = rf = (I/(TT/47rmpC'^ro)'^'*, 
rj,h = rs, and the ratio L°,^T/L^h sa.t urates to eeCd (see 
iMeszaros et al. ll2002l : |Pe'er et al. ||2007| . for further details). 
Higher value of ?7 > 77* does not change this result. 



2.2 



Initial expansion radius, photospheric radius 
and Lorentz factor 



As was shown bv lPe'er et al. I (|2007|), the outflow parame- 
ters, in particular the Lorentz factor rj, the initial expansion 
radius, ro and the photospheric radius, rph can be inferred 
directly from studying the thermal component alone (for 
bursts with known redshift). For completeness of the anal- 
ysis, we bri e fly re peat here the main arguments given by 
iPe'er et oZ.1 (|2007l ). 

The observed temperature of thermal photons emitted 



from the photosphere is T"* ~ 1.4:8riT'{rph)/il+z)E- Due 
Lorentz aberration, the ratio {F^'j^/aT' 



ob4,i/2 



to 



is proportional 



to the photospheric radius divided by the Lorentz factor. 



7^; 



crT"*-" 



1/2 



= e 



jl + zf rph 



di 



V 



(4) 



Here, a is Stefan's constant, ^ is a geometrical factor of order 
unity, d_t is the luminosity distance and f-^^ — L^h/i-rrdj^ 
is the observed thermal flux. Using the parametric depen- 
dence of rph from equation [l] in equation 14] one obtains the 
asymptotic Lorentz factor, 



V 



^{l + zfdL 



2mpC^TZ 



1/4 



L 

T oh 



1/4 



(5) 



Combining this result with the equation for the observed 
thermal flux, L^h ~ (^/2)(7'ph/'"a)~^ , one obtains the ini- 
tial expansion radius, 



ro ■■ 



drTZ 



C(1 + 2)V \2L^t 



-3/2 



and the photospheric radius, 

1/4 



rph 



dlF^iarTZ'' 



^3(l + 2)62mpC3 



r oh 



1/4 



(6) 



(7) 



The values of 77, ro and rph are thus fully determined 
by the observed quantities of the thermal emission, up to 
the uncertainty in the luminosity ratio L/L'fh ^ 1- Further- 
more, equations [5] and [7] imply that the values of 77 and rph 
are not very sensitive to the uncertainty in this ratio. 

Constraining the ratio of the total luminosity released 
in the explosion to the luminosity emitted as thermal pho- 
tons, L/L'fhj can most easily be done if an independent 
measurement of the kinetic energy exists. Such measure- 
ments are provi ded by studying the emission during the af- 
terglow phase (|Wiiers fc Galama 19991: Frail et al. \ 1 2001 : 



iPanaitescu fc Kumar 200ll: iFreedman fc Waxman 2001, 
I Bloom et al. 200S: lBerger et al. 2003; Nysewandcr et aL\ 
[2OO9I), which provides good estimates of the kinetic energy 
remaining after the prompt emission phase. Fortunately, 
such measure ments become ubiqu itous, and are available for 
GRB090902B jCenko et al. 1120111 ): see further discussion in 
^ below. 

Even if afterglow measurements do not exist, the lu- 
minosity ratio L/L'fh can still be constrained indirectly, 
in three independent methods. First, by fitting the non- 
thermal part of the spectrum, one obtains a constraint on 
edee{L/L^h) (see eq. [3]). Since edCe < 1, a lower limit on 
L/L^h is obtained. Second, variability time measurements 
(if they exist) can constrain the initial expansion radius, 
since St ^ ro/c, which, in turn, provides an upper limit on 
the ratio L/L^h via equation [B] And finally, a high value of 
L^h implies that the luminosity ratio should not be high, in 
order to avoid an energy crisis. 



^ Note that in fact at r > rph the kinetic luminosity is slightly 
higher due to energy conversion from the photons above the sat- 
uration radius. While the full treatment will be given below, in 
the content of equation |3] this has little effect, and is omitted for 
clarity. 



3 For photons emitted along the line of sight, the Doppler shift 
is X'(0 = 0) = 2ri. The numerical factor 1.48 results from angular 
integration. 



Connection between thermal and non-thermal emission in GRBs 5 



2.3 Constraint on the emission radius of the 
non-thermal photons 

Observations of high energy ( ~ 10 GeV) photons by Fermi, 
are commonly used in the literature to constrain the Lorentz 
factors of GRB outflows, using the opacity argument. We 
point out though, that the constraints set in the literature 
are often based on an additional assumption, that is that 
the emission radius of the high energy photons, r^ is con- 
nected to the Lorentz factor via r-y = 2T}^cSt, where St is the 
variability time of the inner engine activity. This assump- 
tion, while true in the internal shocks model scenario, has 
two main drawbacks. First, it assumes an a-priori knowl- 
edge of the variability in the Lorentz factor, namely Arj ~ rj 
(this assumption translates into the numerical coefficient); 
and second, it relies on an assumed knowledge of the phys- 
ical variability time, St, which is difficult to be measured 

accurately. 

A different approach was suggested bv IZhang fc Pe'er I 

l|2009f) : by releasing the requirement r-y = 2rfc5t, it was 
shown that the opacity argument can be used to provide 
general constraints in the r.y — rj space. While this method 
does not provide directly the value of 77 (or of r^,), its main 
advantage is that it is not sensitive to the uncertainties men- 
tioned above. Here, w e take the arguments presented by 
IZhang Sz Pe'er I l|2009l 'l one step forward. We first use the 
analysis of the thermal component to estimate the Lorentz 
factor rj. At the second step, we use the constraints found by 
the opacity argument in the r-, — -q plane to deduce a lower 
limit on the emission radius of the non-thermal photons, r^. 
The calculation is performed as follows. The cross sec- 
tion for pair production of photon with energy ei is the high- 
est for interactions with photons of energy £2 ~ {meCp'Y /e\. 
Therefore, considering the Lorentz boosting, a photon ob- 
served with energy e^ax ~ lOeJjfax.io GeV interacts mainly 
with photons having energies 



oh ^ 



{nieC 



2^%2 



26 



.(H-2)2 {1+z) 



2( ob \ 

■ »73(,emax,10j 



MeV. 



(8) 



This energy is about two orders of magnitude higher than 
the energy of the thermal photons (see eq.[2|). We therefore 
do not expect the thermal photons to play a significant role 
in constraining the emission radius of the most energetic 
photons seen by the Fermt-hKT . 

The observed spectrum above a few MeV is often mod- 
eled by a single power law, dN"'' /dt°^dAde°^ = fae°^~" ■ For 
such a spectral fluenc e, the optical d.epth for pair produc- 
tion can be written as (Krolik fc Pier Ill99ll: Woods fc Loeb I 



ll995l : lLithwick fc Sari II2001I : IZhang fc Pe'er Il2009h 



<a> d\ 



/oAt°GcV 



(1 + ^) 



{m^c^ 



(ifl)' 



.(9) 



Here, < cr > is the cross section averaged over all angles; 
for fiat energy spectrum (a = 2), one o btains < a >~ or/ 8 
l|Svensson I I 1987: 'Gupt a fc Zhang llJOOaf l Further note that 
At QoV in equation[9]represents the observed time bin during 
which high energy photons are seen, and thus does not nec- 



* Note that this value is about twice the value presented in 
iLithwick fc Sari I 1I2OOII') . 



essarily correspond directly to the uncertain physical vari- 
ability timejj 

For a flat energy spectrum (a = 2) observed between 
Emin and emajj, the observed fluence is related to the (non 
thermal) luminosity via /o = L^y/47rd| log(£max/emin). 
The requirement that the optical depth to pair production 
of the highest energy photon seen is smaller than unity, is 
translated into a lower limit on the emission radius, 



<<T>L° 



pAt° 



47r log(einax/£min) 



1/2 



3 X 10^^ Lgf Af 



ob 1/2 oh 1/2 -1 

GoV.O Emax.lO '/S '"^"^^ 



(10) 



where log(einax/emin) ~ 20 was taken. For parameters char- 
acterizing GRBs seen by the LAT, the emission radius of 
the non-thermal photons is about 3-4 orders of magnitude 
larger than the photospheric radius (eq. [l]), indicating that 
the observed spectrum must be emitted from (at least) two 
separate regions. 



2.4 The observed non-thermal spectrum: 

additional constraints on the free model 
parameters 

The dissipation at r-y, regardless of its exact nature, pro- 
duces a population of energetic electrons. The energetic elec- 
trons emit the non-thermal part of the spectrum, by radiat- 
ing their energy. There are three main radiative mechanisms 
responsible for the non-thermal emission: synchrotron emis- 
sion, synchrotron-self Compton (SSC) emission and Comp- 
tonization of the thermal photons. 

In order to estimate the relative contributions of the 
different emission mechanisms to the observed spectrum, we 
proceed in the following way. 



2.4-1 Electron energy loss by Comptonization of the 
thermal photons 

The thermal photons serve as seed photons for Compton 
scattering by the energetic electrons. Their existence there- 
fore contributes to the high energy part of the non-thermal 
spectrum. 

The power emitted by Comptonization of the thermal 
photons (in the Thomson regime), relative to the power 
emitted as synchrotron radiation, is given by the ratio of 
the energy density of the thermal photon field to the en- 
ergy density in the magnetic field. At the photosphere, 
the (comoving) energy density of the thermal photons is 
UTh{rph) = aT'{Tj,hf = (L/87rr?C7)2)(rph/r,)-*/^ Since 
above the saturation radius the energy density of the pho- 
ton field drops as uth (xr~^, at the dissipation radius r^ it 
is equal to UTh{r.y) — {L/8Tvr?iri^c)(rph/rs)~'^'^. The energy 
density in the magnetic field assumes a fraction en of the co- 
moving energy density, L/8Tvr?^ri^c. One therefore concludes 



^ E.g., A^Q^y could correspond to the integrated time over sev- 
eral distinct events. 

^ This assumption is taken here as it provides a good first order 
approximation to the observed high energy spectrum. A more 
accurate calculation considering the exact power law index will 
be used in ij3] 
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that the ratio of the energy densities in the thermal photon 
and magnetic field is given by 



Y 



Ub 






-2/3 



3l; 



2/3 8/3 2/3 
I ^3 ''0,8 ' 



(11) 



It is thus clear that for parameters characterizing GRBs, 
the role played by Comptonization of thermal photons is, at 
the least, comparable to the role played by the synchrotron 
emission as a source of energy loss of the energetic electrons. 



2.4-2 Synchrotron spectrum 

The dissipation process is expected to produce a power 
law distribution of energetic electrons with power law in- 
dex p ~ 2.0, above a characteristic Lorentz factor 7™ — 
ee{mp/me) = 184ee,-i.Lj Comparison of the cooling time to 
the dynamical time implies that the entire electron pop- 
ulation is in the fast cooling regime: the cooling time is 
equal to the dynamical time for electrons having Lorentz 
factor^7c = {3TvmeC^ri^r-y)/[aTeBL{l + Y + Y)] = 3.5 (1 + 

Y + Y)~^ L'^^ vi '''1,15 ^B -1- Here, Y has its usual mean- 
ing as the ratio of SSC to synchrotron radiated power, 

Y = PSSC/Psyn. 

The peak of the synchrotron emission. 



-06 



= (3/2)Rr?/(l -I- z){qB-fl/m,c) = 1.5(1 + 



z)-'Ll'\''' 



^54, ^S,-l''7,15^e, 



^ '^ 1 keV, is below the threshold 
energy of the Fermi-GUM detector. Here, q is the electron 
charge and B is the magnetic field at the dissipation radius. 
At the other end of the energy spectrum, comparison 
of the cooling time to the acceleration time provides an 
estimate of the maximum Lorentz factor of the accelerated 
electrons (assuming high efficiency in the acceleration 
process), 7max ~ [67rg/(TT-B(l + Y + Y)]^^^. Synchrotron 
photons emitted by these electrons are expected at energies 
e°r^.^,syn = 240(1 + z)-^ Tja {1 + Y + y)-iGeV, above 
the threshold energy for pair production, and above the 
maximum photon energy seen so far by Fermi. These results 
imply that synchrotron emission is expected to contribute 
to the spectrum at the entire spectral range covered by 
Fermi. 

The fast cooling of the electrons imply that: (I) virtually 
all of the dissipated energy given to the electrons is radiated; 
and (II) for power law index p « 2.0, a flat energy spectrum 
(uFiy oc i^") is expected from synchrotron emission over the 
entire energy range covered by Fermi. Note though that a 
high energy cutoff resulting from pair production can limit 
the maximum observed photon energy to values lower than 

Emax.si/n (sCC ij2]3] aboVc) . 



2.4.3 Comptonization 

There are two sources of Comptonized spectrum: Comp- 
tonization of the thermal photons, and SSC. At low energies, 
below the thermal peak, Comptonization is not a signifi- 
cant source of photons. Hence, the spectrum below the ther- 
mal peak is dominated by synchrotron emission. However, 
above the thermal peak, the three emission mechanisms- 

'' If the dissipation results from shock waves crossing, this equa- 
tion implies a mildly relativistic shock Lorentz factor, Fs — 1 ~ 1, 
and is thus consistent with the internal shocks model. 



synchrotron, SSC and Comptonization of the thermal pho- 
tons contribute in parts to the spectrum. As shown in equa- 
tion [TT] and is further discussed below, the relative contribu- 
tions of the different emission mechanisms are of the same 
order of magnitude (in other words, both Y and Y are of the 
order unity). As a result, it is difficult to determine a single 
dominant emission mechanism at high energies. A conse- 
quence of this, is that the observed spectral index at high 
energies (at the LAT band, above the thermal peak) can- 
not be directly related to a power law index of the energetic 
electrons. 

Comptonization of the thermal emission. Since 
{')mT')/me(? — 0.131/^4 ^ 773' rg^g ee,-i < 1, Comptoniza- 
tion of the thermal emission is in the Thomson limit. The 
Comptonized thermal photon spectrum has characteristic 
breaks similar to the SSC spectrum. Since the electrons are 
in the fast cooling regime, the Comptonized spectrum is 



ob > 



_IC.ob 



expected to rise above 7c r° ~ r° up to a peak at e, 
(4/3)7^(2.87°") « hQ{l^^y2j[^'"\Trl{!el_^GeV, 
roughly as vF^ ex v^^'^ l|Sari fc Esin 1120011 ). At the highest 
energies, photons annihilate by producing pairs, a phe- 
nomenon which can explain the lack of detection of the 
50 GeV photons so far. 

SSC. The ratio of SSC to synchrotron luminosity is 
given by the parameter Y . For e^ ^ tB, the value of Y can 
be approximated asY ^ (te/es)^''^ (|Sari fc Esin 11200 ih FI A 
significant high-energy non-thermal part, as is seen by the 
Fermi-LAT, implies (via eq. [31), high value of te, close to 
equipartition. However, value of es > lO"'^ as is inferred in 
many GRBs (see below), guarantees a value of y of a few 
at most. 

The SSC spectrum rises as vFi, oc v^'^ below the peak 
of the SSC emission, which is expected at el^'"^ = 27^5^ « 
100(1 + z)-^L\{^e]l'^_^et_-^r-\r^MeV. At higher energies, 
the spectral shape follows a similar power law as the syn- 
chrotron spectrum, i.e., a flat energy spectrum is expected 
for p ~ 2. 

The rise parts of both the thermal Comptonization and 
the SSC emission are independent of the power law index 
of the accelerated electrons. For power law index p ~ 2.0 a 
flat, or slightly decaying synchrotron (energy) spectrum is 
expected at the entire Fermi energy range (see H2.4.2l above'). 
The combined effects of the flat synchrotron emission with 
the rise of the Comptonized spectrum therefore results in a 
mild increase in the total observed spectrum at high ener- 
gies, above the thermal peak (see ij3] below). Demonstration 
of spectral decomposition into its basic physical ingredients 
is presented in ij^ 

The rising of the SSC component below ~ 100 MeV, 
combined with the fact that the synchrotron spectrum is 
expected to be flat (or slightly inverted) and that Y is not 
expected to be much larger than a few, imply that at low 
energies (below the thermal peak), Comptonization is not 
expected to play a signiflcant role. The main emission mech- 
anism below the thermal peak is therefore synchrotron emis- 
sion. Since the Fermi-GBM detection range is above ej^, the 
spectrum below the thermal peak is expected to be sensi- 
tive to the power law index of the accelerated electrons, p. 



Note that HY = Y, one obtains Y ~ (es/2ee)^/^. 
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Thus, measurement of the flux at low energies can provide 
an indication for both the values of es and of p. 

Finally, we point out that the analysis carried in this 
section holds only as long as r-y ^ r^h- For r^ ~ r^h, the 
effect of Comptonization is much more complicated due to 
the fact that the electrons cooling time is oc r^. Therefore, 
for small dissipation radius, the cooling time is much faster 
than the dynamical time for all electron energies, and 7c ~ 1. 
For such rapid cooling, additional physical phenomena, that 
are not considered here, become important. Direct Compton 
scattering of the energetic photons provide the main source 
of heating, resulting in a q uasi steady state d istribution of 
mildly-relativistic electrons IjPe'er et al. II2005I ) . Close to the 
photosphere, multiple Compton scattering by electrons in 
this quasi steady state produces a flat energy spectr um for 
a large region of parameter space (jPe'er et al. II2006I ). 



3 DEMONSTRATION OF THE ANALYSIS 
METHOD: GRB090902B AS A CONCRETE 
EXAMPLE 

The bright, long GRB090902B, w hich is one of the bri ghtest 
bursts observed by LAT to date (jAbdo et al. Il2009al ). pro- 
vides an excellent example for demonstrating our analysis 
method. This is because of two very pronounced proper- 
ties of its prompt emission spectrum. First, time resolved 
spectral analysis reveals a signiflcant power law component 
in the LAT data (emission was observed up to 30 GeV ), 
w hich is clearly distinct from the usual "Band" function used 
bv lAbdo et al. I (12009^1 t o fit the spectrum in the sub-MeV 
range (JRvde et al. ll2O10l : IZhang etal. Il201lh . Moreover, the 
fact that the high energy photons were delayed with respect 
to the photons at the sub MeV peak indicates a different 
origin. Second, the "Band" function used in fltting the sub 
MeV peak is extremely steep on both sides (low energy spec- 
tral slope nie) oc e" , with a « — 0.5..0.3, and high energy 
spectral slope /3 ~ — 3.. — 5), resulting in an unusually nar- 
row peak. The steep spectral slopes seen in the sub-MeV 
range make it impossible to flt this spectrum with a model 
that contains only synchrotron and SSC. 

On the other hand, both properties of the prompt spec- 
trum fit perfectly into the framework suggested here: first, 
the sub-MeV peak is naturally modeled with the thermal 
emission component. While a single black body provides a 
sufficient fit, and is used to deduce the values of ri,ro and 
Tph, better flts are obtained with multi-color black body, as 
expected from a theoretical point of view (jPe'er fc Rvde I 
[201l|). Second, the non-thermal part (the high energy power 
law which extends to lower energies), can be easily explained 
by a combination of synchrotron, SSC and Comptonization 
of the thermal photons. Moreover, by doing so, we deduce 
the physical properties in the emission sites of both the ther- 
mal and non-thermal components, hence we obtain a com- 
prehensive physical picture of the properties of t his burst. 

In the original analysis, lAbdo et al. I (|2009al ) separated 
the observed prompt emission into several time bins. The 
most restrictive constraints on the emission radius of the 
non-thermal photons arise in time interval (c), 9.6 — 13.0 s 
from the trigger, since in this time interval an 11.16 GeV 
photon was observed. We therefore focus our analysis on 
this time interval. 



3.1 Analysis of the thermal component: Lorentz 
factor, photospheric radius and dissipation 
radius 

In order to deduce the value of the Lorentz factor by using 
the method presented in i)2.2l one needs to flt the sub-MeV 
peak with a single black body spectrum, and study its prop- 
erties (temperature and thermal flux). While this can be 
done for the data in time interval (c), since the properties 
of the outflow vary on a shorter time scale (variability time 

of ~ 0.1 s was observed), smearing of the black body spec- 
trum is expected. The shorter the time interval used for the 
fits, the higher the quality of the blac k body fits obtained 
(|Rvde et al. II2OI0I : [Zhang ef al. Il2011^ . On the other hand, 
shorter time intervals result in lower quality of the high en- 
ergy data, which is sparse. 

We thus use time interval as short as possible to flt the 
narrow peak with a single black body, in order to deduce 
the hydrodynamics of the flow, and in particular obtain the 
Lorentz factor. We then use the longer time interval (the 
full time interval (c)) to study the properties of the non- 
thermal part of the spectrum. For the long time interval, we 
flt the spectral peak with a multic olor black body, which, 
as was shown bv lRvde et al. I l|2010f ), provides better flts to 
the peak. Clearly, by choosing to fit the longer time inter- 
val during which the values of the parameters vary, we lose 
an accuracy in the flts. However, we stress here, that our 
purpose in this paper is not to provide the best statistical 
flts (in terms of x^) to the data. Instead, our goal here is to 
prove that one can obtain an acceptable flts to the data in 
the entire Fermi energy band, in the sense that the curves 
plotted all fall within the ilcr error bars of the data points. 
By doing so, we show that we are able to provide an ac- 
ceptable, complete physical interpretation to the observed 
data, although, necessarily, we are not able to capture many 
second-order effects. 

For t he shorter time int erval we use the results pre- 
sented by iRvde et al. I ((20101), who showed that the nar- 
row time interval, 11.008 — 11.392 s from the trigger, dur- 
ing which the 11.16 GeV photon is seen, provides suffl- 
cient data to analyze the sub-MeV spectral peak. Dur- 
ing this time interval, the sub-MeV peak can be flt- 
ted with a single black body, with observed tempera- 
ture T°^ = 168 keV, and observed thermal flux F^i = 
1.96 X 10^^ ergcm^^ s~^. Equation 0] thus implies a ra- 
tio 7^ = 1.55 X 10 "^^ cm. At redshift z = 1.822 
ICucchiara et al. II2009I ). the (isotropic-equivalent) thermal 
luminosity is L^^ = AivdlF^l ^ 4.6 x lO^^ergs" ^ Using 
^ = 1.06 in equations [5l[6l and 171 (|Pe'er et al. |[2007h . one ob- 
tains 7? = 764(L/L??;,)l''^ rph = 6.0 x lO" {L/L^\y^^cm 
and r-o = 7.9 x 10** {L/2L^\)-^^^ cm. 

Estimating the ratio {L/L^j^) is done in the following 
way. Due to the rapid cooling of the electrons, nearly 100% 
of the energy that is converted to the energetic electrons 
in the dissipation process is radiated in the form of non- 
thermal photons, that is L?^t — CeedLk. Here, Lk is the 
energy available in kinetic form above the photosphere. By 
fltting the non-thermal part of the spectrum, one obtains 
L°j}t ^ 0.9LTh, or ee^d - 0.9(Lrh/Lfe) (see Figure[lJ. These 
flts are best done using a numerical simulation, since part 
of the dissipated energy is released outside the observed en- 
ergy range of Fermi. After the main dissipation, the remain- 
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ing kinetic luminosity, Lag ~ Lk{l — td) is the available 
lu minosity for the afte rglow emission phase. As was found 
bv lCenko et al. I (J201lh . the energy release during the after- 
glow phase is ~ five times less than the energy release during 
the prompt emission phase. However, since the luminosity in 
time interval (c) is about twice the average luminosity dur- 
ing the prompt phase, we estimate the luminosity (thermal 
-I- non-thermal) in time interval (c) to be about ten times 
higher than the luminosity during the afterglow phase. 



Lag 



T oh I T ob 



Lkjl - (-d) 

L°^\+Lke,ed 



10' 



Using eettj ~ 0.9(Ly^/Lfe), one obtains the relation 
9 



Ed = 



1.9ee + 9 



(12) 



(13) 



This result immediately implies very high dissipation effi- 
ciency, since for any value of ee, e^ ~ 0.83. 

At the photosphere, the thermal luminosity is equal 
to L^Ti = {L/2)(r-ph/rs)~^ '^ , and the kinetic luminosity is 
equal to Lk = (^/2)[2 - (rph/r^)"^''^]. At larger radii, the 
kinetic luminosity is assumed unchanged, up until the dis- 
sipation radius. Using the observed relation L^^/L^ti = 
€dteLk/L^^ ~ 0.9 and the results obtained in eauation ll3l 
one obtains the luminosity ratio , 



(1.19ee+0.9) 



r ob 



(14) 



For equipartition value, e^ — 0.33, one obtains L/L^ti = 
3.9, which imply r) — 1070 (and e^ ~ 0.94). However, 
we consider this value of the Lorentz factor as an upper 
limit, due to two reasons: first, this result implies very high 
total GRB luminosity, L = 1.8 x 10''^ ergs~^; and sec- 
ond, using equation El it implies an initial expansion ra- 
dius ro — 3.0 X 10* cm, which translates into very short 
variability time, St — vq/c = 10 ms. Lower value of e^ re- 
sults in higher Lorentz factor, but also leads to significantly 
higher total luminosity, significantly lower variability time 
and nearly 100% dissipation efficiency, which we consider 
unlikely. 

Using the extreme value e^ — 1, one obtains L/L^f^ = 
2.09, ri = 920, e^ = 0.83 and ro ~ 7.4 x 10® cm, which 
translates into physical variability time St ~ 25 ms. We thus 
conclude, that the asymptotic value of the Lorentz factor is 
in the range 920 ^ ry ^ 1070, and the value of e^ is at 
or slightly above equipartition. We further deduce that the 
photospheric radius is rph — 7.2 — 8.4 x 10^^ cm, and that 
the dissipation efficiency is e^ « 85% — 95%. Using these 
values of rj in equation 1101 leads to the conclusion that the 
emission radius of the non-thermal photons is constrained, 
r-., ^ 3.5-4.1 X lO^-'^cm. 

The derived value of ro — 3.0 — 7.5 x 10* cm implies 
ro ~ 33 — 80 rjsco, 10, where r/sco, lo = QGM/c^ is the 
inner most stable circular orbit (ISCO) of a non-rotating 10- 
solar mass black hole. As the exact mass of the progenitor 
of GRB090902B is unknown, we can deduce that ro is of 
the order of few tens of the fSCO radius. According to the 
theory adopted here, ro marks the initial expansion radius. 
As the theory of jet acceleration is not fully developed yet, 
the values obtained may thus be used to constrain models of 
jet acceleration to this scale. Alternatively, we note that this 
value is very close to the value obtained in several numerical 



models (e.g.. lAlov et al. II2OO0I : IZhang et al. II2003I '). While 
in these models the jets are assumed to be produced closer 
to the ISCO radius, the acceleration begins only at larger 
radii due to the fact that at smaller radii the jet is not well 
coUimated. 



3.2 Numerical calculations: further constraints on 
the free model parameters 



As discussed in i]2.4l above, following the dissipation pro- 
cess, simple analytical descriptions are insufficient to de- 
scribe the spectrum, due to the roughly similar contribu- 
tions from the different physical phenomena (synchrotron, 
SSC and thermal Comptonization), as well as the cutoff 
resulting from pair production. Therefore, in order to de- 
rive the spectral dependence on the different values of the 
free parameters, as well as confirm the analytical calcula- 
tions presented above, we calculate numerically the pho- 
ton and particle energy distribution for the different sce- 
narios. In our calculations, we use the time-d e pende nt nu- 
merical code presented in IPe'er fc Wajcman I (|2005h . This 
code solved self-consistently the kinetic equations that de- 
termine the temporal evolution of e* and photons, describ- 
ing cyclo-synchrotron emission, synchrotron self absorption, 
direct and inverse Compton scattering, pair production and 
annihilation and the evolution of high energy cascade. 

This code has two great advantages, which make it ideal 
in the study of the prompt spectra. First, it has a unique 
integrator, that enables solving the rate equations that gov- 
ern the time evolution of the particles and photons energy 
distribution over the entire energy range. The code is able 
to calculate simultaneously processes happening over more 
than 15 orders of magnitude in time and energy ranges, 
thereby covering the entire spectral range, from radio up 
to the TeV band. The second advantage of the code is a 
full treatment of the various physical process, including the 
full cross sections (e.g., Klein Nishina effect is inherently 
taken into account; or that for mildly relativistic electrons, 
the full cyclo-synchrotron emission spectrum is calculated; 
etc.). 

In modeling the spectrum in time interval (c), the 
sub MeV peak is best described as a mult icolor black 
body spectrum, with varying amplitudes l|Rvde et al. I 
EoiO): F^l(v-T^^^) = ]^2^^{dA/dT)B,{T)dT, where 

B^{T) = [2h/c^)v^/[e'"'/'^ ~ 1) is Planck function. 
The amplitude A{T) is temperature dependent, A{T) = 
^(rmax)(rmax/r)^~'', and is normalized such that the to- 
tal flux is equal to the observed flux at the spectral peak, 
Fyjj — 1.82 X 10~* ergcm"^ s~^. The normalization con- 
stant q and Tmax are determined by fitting the spectrum, 
Tmax = 328 keV and q — 1.49. Tmin cannot be determined, 
and its exact value is unimportant for the fits, as long as 

-^ mill ^^ -^niax. 

While the best value for q is found by fitting the data, 
we note that this multi-color description of the thermal part 
of the spectrum is inherent to emission from relativistically 
expanding plasma. At any given instance, an observer sees 
simult aneously pho tons emitted from a range of radii and 
angles l|Pe'er Il2008l 'l. The Doppler boosting of photons emit- 
ted at high angles to the line of sight is smaller than that 
of photons emitted on the line of sight. Therefore, a pure 
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Planck function in the comoving frame is inevitably observed 
as multi-color black body. 

In the framework developed here, the value of q is re- 
lated to the spectral index via F^ oc v"^'^ . This can be 
seen by noting that dA/dT oc T''"^, and hence Fth{v) oc 
J dT{dA/dT)B„{T) oc u''~^; the last equality is easily ob- 
tained by replacing the integrand from T to z — hv/T. Thus, 
the value of q found by fitting the data implies spectral index 
F^ oc i/"*^. This index is significantly softer than the index 
expected for pure black body (Fi , (X v^ ), or from the sce- 
nario considered bv lBeloborodov I ((20101), of fixed comoving 
temperature, in which F^, oc v^'^. On the other hand, this in- 
dex is harder than th e expected index at late times, F„ oc i/" 
l|Pe'er fc Rvdel201ll 'l. At these times, off-axis emission dom- 
inates the spectra. When the spectrum is dominated by on- 
axis emission, in the spherically symmetric scenario it is ex- 
pected to be Fi, oc u^ [Lundman et. al., 2011, in prep.]. The 
fitted spectral index is thus in between these values. We 
find this result encouraging, given that (I) here the fitted 
spectrum is integrated over a finite time interval, hence we 
average over spectra obtained at different times, and (II) the 
theories are developed for the 'pure' case of constant out- 
flow Lorentz factor. Hence, the results obtained are in good 
agreement with the theoretical expectations|j 

In producing the spectrum, we assume that at radius 
r^, a fraction ed of the kinetic energy is being dissipated (by 
an unspecified dissipation process). The energetic electrons, 
which assume a power law distribution with power law index 
p, carry a fraction te of the dissipated energy, and the mag- 
netic field carries a fraction tB of this energy. The multicolor 
black body spectra serve as background spectra for all the 
various processes (mainly Compton scattering by the ener- 
getic electrons, but also other processes such as, e.g., pair 
production). The code tracks the evolution of the spectrum, 
during the dynamical time. Here, we present the results at 
the end of time interval (c), i.e., we assume that the pro- 
cesses take place during an observed time of 3.4 seconds. 



3.3 Numerical results 

The numerical fits to the spectrum of GRB09092B at time 
interval (c), 9.6-13.0 seconds after the trigger, are presented 



in Figures \T\ - l4l^°l In Figure [TJ we demonstrate the lin- 
ear dependence of the non-thermal flux on the value of eg. 
For the fits shown in this figure, we chose parameter val- 
ues that fulfill the requirements in the previous sections 
and are typical for GRBs. Thus, we chose large dissipa- 
tion radius, Tj = 10^^ cm, strong magnetic field, es = 0.1, 
and electron power law index p — 2.0. The Lorentz factor 
was chosen to be F = 910, which implies L/L^f^ = 2.0, or 
L = 9.2 X 10^^ ergs~^. These values imply, via equation 1 141 
ee = 0.5 and ea = 0.9. The fit to the Fermi (GBM -I- LAT) 
data at time interval (c) appears as the solid (blue) line 
in Figure [l] The fit shows the combines spectrum resulting 



^ Addition al broadening may result from sub-photospheric dissi- 
pation, see lRvde et al. I 1120111 ) 

^^ These can be compar ed to the spectra presented in Figure 
3 in lAbdo et al. I 1120093^ , although note that in this work the 
spectrum is presented at a different time interval. However, the 
spectral shape and the main spectral features (such as the ratio 
of the flux at the peak to the flux at ~ lOkeV) are similar. 



from both the thermal peak (at ~ MeV), and the nearly 
flat energy spectrum [lyFi, oc i^") resulting from synchrotron 
emission from electrons in the fast cooling regime. We fur- 
ther added a scenario in which e^ is three times smaller 
(te = 0.17, dashed green line), which demonstrates the lin- 
ear dependence of the non-thermal flux on the value of e^. 
We point out that although a power law index p = 2.0 was 
chosen, the combined effects of flat {vF^, oc u°) synchrotron 
spectrum and rising [lyF^, oc i/^") Comptonization spec- 
trum, lead to a slight increase in the high energy spectral 
slope, which is consistent with the slope seen with Fermi. 

We further note that the fit to the Wien part of the 
thermal component falls slightly below the ±lcr error bars 
of the data (the shaded, yellow areas in the figures). This 
discrepancy can be easily understood as due to smearing of 
the data: the data presented in the figures is averaged over 
several seconds, during which the properties of the outflow 
(such as the Lorentz factor) slightly vary, while in the numer- 
ical flt we assume steady values of the physical parameters. 
Variation in the parameters values inevitably lead to smear- 
ing of the signal, which is translated to a high energy decay 
which is somewhat shallower than the exponential cutoff of 
the thermal spectrum considered by the fits. Indeed, detailed 
ana lyses of time resolved spectra done bv lRvde et al. I (|2010l ) 
and lZhang et al. I (|201ll ) show that as the time interval con- 
sidered becomes shorter, the exponential decay above the 
thermal peak becomes more and more pronounced, hence 
the (multi color) black body function used in fitting the 
peak becomes better the shorter the time bin is. Nonethe- 
less, as explained above, we chose here to fit the data in the 
entire time interval (c), since reducing the time interval re- 
sults in poor quality of the high energy part of the data (the 
non-thermal part). We therefore find it appropriate to use 
the fits presented, as they serve the main goal of this paper: 
to demonstrate that the physically motivated, hybrid (ther- 
mal -I- non-thermal) model provides acceptable fits, which, 
moreover, enable a good estimate of the physical conditions 
at both emission sites. 

In Figure[21 we consider different dissipation radii: r^, — 
10^^ cm (solid, blue), 10^*^ cm (dashed, green), 10^^'^ cm 
(dash-dotted, red) and 10^^ cm (dotted, purple). As the nu- 
merical code considers the full cross section for pair produc- 
tion, the numerical results are more accurate than the an- 
alytical approximations presented in i]2.3l and can be used 
to validate them. The results presented in Figure [2] indeed 
confirm the main conclusion obtained analytically, that is 
that the observation of the 11.16 GeV photon necessitates 
the dissipation radius to be above 10^^'^ cm. 

At larger radii, the high energy non-thermal part of the 
spectrum is not very sensitive to the exact dissipation radius. 
As shown in Figure [5J for dissipation radii r^ ^ 10^^'^ cm 
it is possible to obtain numerical results which are within 
±1(7 errors of the empirical "Band" fit. In order to achieve 
this, high value of es and a slight tuning of the value of Ss 
is required. Thus, for r^, = 10^® cm, a value of ee = 0.4 was 
chosen, while for the other fits in this figure, ee = 0.5 was 
found adequate. 

The numerical results show that the dependence of the 
pair production cutoff at high energies on r^, is weaker than 
the analytical approximation presented in equation 1101 The 
main reason for this is the contribution from the thermal 
photons, which is neglected in the derivation of equation 
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Figure 1. The dependence of the non-thermal flux on the fraction 
of energy given to the electrons, te. The data of GRB090902B arc 
from the Nal (0,1), BGO(0,1) and LAT (back and front) detectors 
at time interval (c), 9.6 - 13.0 seconds after the GBM trigger. 
The light yellow shaded ar ea show the ilcrfit to the data at 
this time interval, ma de bvlAbdo et al. I ll2009al V Note that this 
area is calculated by lAbdo et al. I l|2009a|) by fitting a "Band" 
plus a single power law spectra, hence the apparent discrepancy 
between the high energy (LAT) data points and the shaded area. 
In our work, the parameter values chosen are: dissipation radius 
r-Y = lO^'^ cm, bulk motion Lorentz factor F = 910, power law 
index of the accelerated electrons p = 2.0, ORB luminosity L = 
9.2 X 10®^ ergs~^ and fraction of dissipated kinetic energy e^ = 
0.9. Shown are the simulation results for e^ = 0.5 (blue, solid 
line), and te = 0.17 (dashed, green line). The ~ MeV peak is 
fitted with a (multi-color) black body spectrum; the non-thermal 
flux is linear in the value of eg. 



1101 We can therefore conclude, that further constraints on 
the dissipation radius can not be obtained directly from the 
prompt spectrum, without additional assumptions. Finally, 
we point out that the small bump obtained in the scenario of 
r-y — 10^^ cm at TmeC^ /{1 + z) ~ 150 MeV, results from pair 
annihilation process, which is more pronounced at small dis- 
sipation radii due to the more rapid production of e pairs. 
Such a small bump is difficult to be observed, and indeed has 
not been observed so far. Moreover, since the Lorentz factor 
of the flow likely varies during the time interval considered 
here, it is expected to be smeared. 

In Figures [3] and U we examine the dependence of the 
spectra on the uncertain values of the power law index of 
the energetic electrons and the fraction of dissipated en- 
ergy carried by the magnetic field. In Figure [31 we consider 
three values of the power law index: p — 2.0 (solid, blue), 
p — 2.2 (dashed, green) and p = 2.5 (dash-dotted, red). 
As explained in section [2.4.31 above, and is further demon- 
strated in appendix 1X1 the high energy part of the spectrum 
(above the thermal peak) is governed by Compton scattering 
which is in the rising part of the spectrum (below £„ ), and 
therefore the spectrum is not sensitive to the exact power 
law index of the accelerated electrons. We thus conclude 
that for GRB090902B-type bursts, for which the thermal 
peak is pronounced, j j observations at high energies cannot 
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Figure 2. The dependence of the non-thermal flux on the dissi- 
pation radius, r^. We show the numerical results for dissipation 
occurring at r-,, = 10^^ cm (solid, blue), 10^^ cm (dashed, green), 
10^^'^ cm (dash-dotted, red) and 10^^ cm (dotted, purple), su- 
perimposed on the Fermi data and the ±1(T "Band" function fit 
to the data (light yellow shaded area). Values of eg = 0.33 and 
ee = 0.5 were chosen, apart from the fit for rj = 10^^ cm, where 
eg = 0.4 is chosen. All the other parameters are the same as in 
Figure 1. Below 10^^'^ cm, pair production limits the maximum 
observed energy of photons to ~ GeV, and is thus inconsistent 
with the LAT observation of 11.16 GeV photon at this time in- 
terval. 



be used to constrain the power law index of the accelerated 
electrons. On the other hand, the low energy part of the 
spectrum (below the thermal peak) is dominated by syn- 
chrotron emission, and is thus sensitive to the power law 
index of the electrons. Unfortunately, most of the effect is 
expected below the threshold energy of the Fermi- GBM 
detector, and therefore only weak observational constraints 
exist. We can therefore conclude that a power law index in 
the range 2.0 ^ p ^ 2.2 is consistent with the data, and 
even the higher value of p = 2.5 can be consistent, however 
for such high power law index a somewhat fine tuning of 
the other model parameters (in particular, very high value 
of e_B, close to equipartition) is required. 

Examination of the spectral dependence on the value 
of e_B is presented in Figure |31 The three fits presented in 
this figure, equipartition (es ~ 0.33; solid, blue), eg ~ 0.1 
(dashed, green) and eg = 0.01 (dash-dotted, red) show that 
the value of €b cannot be too far below equipartition: low 
value of es results in a too low flux at low energies (below the 
thermal peak), which is inconsistent with the observation. 
It also leads to a more pronounced Compton peak, which is 
seen at the high energies. We can thus conclude that value of 
e_B close to equipartition is needed to be consistent with the 
observations. This high value justifies the need for numerical 
analysis, since although e^ is very high, both Y and Y are 
close to unity. 



Some other bursts, such as, e.g., GRB080916C show much less 



pronounced thermal peak, hence Comptonization may play a sub- 
dominant role for GRB080916c-type bursts. 



Connection between thermal and non-thermal emission in GRBs 11 




10 10 10 

Obs. energy [eV] 

Figure 3. The dependence of the non-thermal flux on the power 
law index of the accelerated electrons. On top of the Fermi data, 
shown are the numerical results for p = 2.0 (solid, blue), p = 2.2 
(dashed, green) andp = 2.5 (dash-dotted, red). Dissipation radius 
r^ = 10^^ cm, eg = 0.5, eg = 0.33 and all other parameter values 
same as in Figure 1 are chosen. The high energy spectrum is 
nearly insensitive to the exact value of p in the range considered, 
2.0 — 2.5. However, the low energy part (below the thermal peak) 
may provide indication for 2.0 ^ p ^ 2.2. 
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Figure 4. The dependence of the non-thermal flux on the frac- 
tion of energy carried by the magnetic field, es- Shown are the nu- 
merical results for equipartition {tg = 0.33; solid, blue), tg = 0.1 
(dashed, green) and eg = 0.01 (dash-dotted, red). Dissipation ra- 
dius r^ = 10^® cm, ee = 0.4 and all other parameter values same 
as in figure 1 are chosen. While the effect of eg on the high energy 
spectrum is minor, the low energy flux (below the thermal peak) 
necessitates high value of eg , close to equipartition. 



4 SUMMARY AND DISCUSSION 

In this paper, we considered the effect of thermal emission 
on the observed GRB prompt emission spectrum. Being a 
natural outcome of the GRB fireball model, thermal emis- 
sion is an inherent part of the prompt emission seen. As we 
showed in i]2.1l (ea.[3)). it is expected to be more pronounced 
for bursts with higher Lorentz factor. As we show here, the 
inclusion of the thermal emission in the calculation of the 
prompt emission spectra, enables one to obtain a complete, 



self consistent physical model of the prompt emission spec- 
trum seen over the entire FERMI energy range. This is in 
contrast to the "Band" function fits, which do not carry 
any physical interpretation, and, in addition, require extra 
models to be able to fit the spectrum at high energies. 

According to our model, the sub-MeV peak often seen 
is interpreted as being composed of multi-color black body 
emission from the photosphere. The non-thermal, high en- 
ergy part seen in several bursts by the LAT instrument 
is interpreted as resulting from combined emission of syn- 
chrotron, SSC and Comptonization of the thermal photons, 
following an episode of energy dissipation that occurs at 
large radius above the photosphere, r-y > rph- Observations 
of high energy photons can be used to constrain the dis- 
sipation radius (eq. llUp . which is clearly above the pho- 
tospheric radius (eq. [!}. Moreover, as we showed in i)2.4l 
the relative contributions of synchrotron emission, SSC and 
Comptonization of the thermal photons, denoted by the pa- 
rameters Y and Y , are expected to be of the same order 
of magnitude. Hence, the three emission mechanisms have 
roughly similar contributions to the high energy (above the 
thermal peak) part of the spectrum. This fact makes it dif- 
ficult to directly determine the power law index of the ac- 
celerated electrons from measurements of the high energy 
spectral slope. 

The separation made here between thermal and non- 
thermal emission, makes it possible to deduce the values of 
the free model parameters. First, by analyzing the photo- 
spheric part, one can deduce the value of the Lorentz fac- 
tor, the initial expansion radius and the photospheric radius 
(see i]2.2p . Then, by analyzing the non-thermal part, one can 
constrain the dissipation radius, r^ (eq. llOfl , and place con- 
straints on the power law index of the accelerated electrons, 
p, and the strength of the magnetic field, eg. By combining 
the fluxes of the thermal and non-thermal parts, one can fur- 
ther constrain the combined fractions of dissipated kinetic 
energy (e^) that is received by the energetic electrons {ce). 
Further separation of the values of these variables is more 
tricky, but can be done with the help of afterglow observa- 
tions [see 3331 eas.[T2l-[T4)). 

We demonstrated our analysis method on the bright, 
long GRB090902B. This burst is ideal for our demonstra- 
tion purposes, because of the clear separation between the 
thermal and non thermal components seen, and the very 
pronounced thermal peak: both above and below the sub- 
MeV peak seen in this burst, the spectral slopes are too 
steep to enable fitting the spectrum with any combination 
of synchrotron and SSC emission models. However, as we 
showed in i|3l Figures [l] - |4l excellent fits are obtained us- 
ing the hybrid (thermal -I- non-thermal) model considered 
in this paper. While we stress again that we did not make 
any attempt to obtain the statistical best fits to the data, we 
are clearly able to obtain fits that are within the ±lcr error 
bars of the data, over a very broad band - about 6 orders 
of magnitude spectral range. These fits are obtained using 
well understood emission mechanisms. We can therefore con- 
clude that, at least for this burst, the 'Band' function is well 
represented by a combination of physical emission processes. 
By doing so, we gain an insight into the physical conditions 
in the emitting regions. 

According to our interpretation, the very pronounced 
~ MeV peak represent the contribution of the thermal com- 
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ponent, and the broad band high and low energy spec- 
trum result from non-thermal processes. By making this 
separation, we were able to constrain the hydrodynamic 
properties of the outflow and the physical properties at 
the emission site: the initial expansion radius was found 
to be To = 2.9 — 7.5 x 10* cm, the photospheric radius 
Tph — 7.2 — 8.4 X 10^^ cm and the Lorentz factor is in the 
range 920 ^ »? ^ 1070. The main source of uncertainty in 
these values is the unknown kinetic luminosity. 

Fit results from the interval during which an 11.16 GeV 
photon was observed, constrain the dissipation radius to be 
r-y ^ 3.5 — 4.1 X 10^^ cm. Combined measurements of the 
thermal and non-thermal parts imply very high dissipation 
efficiency, ed — 85% — 95%, and very high fraction of dissi- 
pated energy carried by the energetic electrons: e^ is at or 
above equipartition value. Measurements of the low energy 
part of the spectrum (below the thermal peak) imply high 
magnetic field, es close to equipartition. The power law in- 
dex of the electrons is more difficult to constrain, as the high 
energy spectral slope results from a combination of flat (or 
slowly decaying) synchrotron part, and rising Comptonized 
part. As such, the observed spectral slope does not directly 
corresponds to the electrons power law index. From the nu- 
merical fits, we concluded that 2.0 ^ p ^ 2.2 are consistent 
with the data, and that even p = 2.5 is marginally consistent 
with the data. 

The numerical results presented in Figure [2] show high 
energy cutoff due to pair production phenomenon. While 
this cutoff is consistent with the analytical prediction in 
equation 1101 it also shows the limitation of the analytical 
approximations, which are commonly in use. Due to the in- 
clusion of the thermal photons, the dependence of the cutoff 
energy on the emission radius is weaker than the simple ap- 
proximation given in equation 1101 

Additional constraints may be obtained by adding addi- 
tional information, albeit with a higher level of uncertainty. 
For example, assuming the internal shocks scenario, and con- 
sidering that the 11.16 GeV photon was seen after a delay of 
11 s from the GBM trigger, implies, for constant Lorentz fac- 
tor and IS.!) ~ 77, dissipation radius of r-y ~ rj^ct « 2 x 10^^ 
cm. We point out that this assumption is consistent with 
the observations, since the much more rapid variability time 
(^ 0.1 s) can be attributed to emission from the photo- 
sphere, rather than the high-energy non-thermal part. 

While GRB090902B shows very pronounced peak and 
clearly separated high energy component, which make it 
ideal for demonstrating our analysis method, this is not 
the case in many GRBs (see, e.g., recent analysis by 
IZhang et ffl"ni201lh . In fact, in many GRBs, the sub-MeV 
peak often seen is not as pronounced as in GRB090902B: the 
low and high energy spectral slopes are not as steep, and so 
in many cases pure single black body spectrum is too narrow 
to fit the sub MeV p eak. In some cases, e.g., GRB100724B 
l|Guiriec et al. 1 120111 ). a weak thermal component, that is 
not associated with the main peak of the spectrum, can be 
identified. In many other cases, there is no clear evidence for 
the existence of thermal component, a s black body spectrum 
is not clearly identified (JNava et al. 112011.1^1 . Moreover, of- 



^^ Although the steepness of the low energy spectral slope can 
be viewed as an indirect evidence. 



ten the high energy spectral slope (above the thermal peak) 
decays with a power law index much shallower than seen in 
GRB090902B. Thus, to date, the spectrum of GRB090902B 
is unique by having such a pronounced thermal peak. 

Clearly, a lack of very pronounced thermal peak in 
most bursts is a major drawback to the ideas raised here. 
There are several ways that can explain these observations. 
One possibility is suppression of the photospheric compo- 
nent which is expected in Poynting flux dominated flow 
(|Zhang fc Pe'er l2009l : rZhang fc Yanl201ll '). In this scenario, 
the magnetic field, rather than the photon field, serves as 
an energy reservoir. As the magnetic energy is gradually 
dissipated, most of the emission occurs at large radii above 
the photosphere, leaving only a weak photospheric signal. 
Within the fireball model itself, the pronunciation of the 
thermal peak strongly depends on the value of the Lorentz 
factor, r, which is a free parameter of the model (see equa- 
tions [51E}. While there is a clear indication for F ~ 10^ in 
GRB090902B and several other burs ts, in many bursts th e 
value of r is lower, ~ 10^'^ (see, e.g.. iRacusin et al. II2OIII ). 
Thus, for these bursts, the thermal emission, while expected 
to exist, is not as pronounced as in GRB090902B. 

Alternatively, these observations can be explained in 
a framework similar to the one used here, namely that the 
spectrum may be dominated by Comptonization of the ther- 
mal photons, following energy dissip ation tha t occu r s clos e 
to the photosphere. As shown bv iPe'er et al. I l|2005l . l2006l '). 
in this scenario, nrultiple Compton scattering by electrons 
at a quasi steady state produces a flat energy spectra for 
a large parameter space region. Since the dissipation radius 
(r-y) can in principle take any value above the saturation ra- 
dius, different dissi pation radii can le ad to very different ob- 
served spectra fsee lPe'er et al. II2OO6I): in particular, high en - 
ergy power law tail is obtained ( Lazzati fc Begelman II2OI0I ). 
Thus, su b-photospheric diss ipation can reproduce the Band 
function (JRvde et al. 1120111 ). In fact, it could very well be 
that the uniqueness of GRB090902B originate from a very 
large dissipation radius, r-, 2> rph. Only under this condi- 
tion one is able to make such a clear separation between the 
thermal and non-thermal parts of the spectrum, which are 
otherwise coupled. 

This possibility can also explain the lack of GeV emis- 
sion in many bursts. If r^ is not much larger than rph, then, 
by definition, the optical depth to scattering is high. Since 
the cross section to pair production is similar to or, the 
optical depth to pair production is high too (see equation 
llOf) . As a result, GeV emission is attenuated. Alternatively, 
attenuation of GeV emission is expected if the outflow is 
highly magnetized: in this scenario, synchrotron emission 
dominates over IC scattering, leading to attenuation of the 
high energy emission. 

In addition to the complex relations between the ther- 
mal and non-thermal parts of the spectra, there are two 
effects which are often being neglected. First, the fits are of- 
ten made to time integrated spectra. As the properties of the 
outflow, in particular the Lorentz factor, vary on a very short 
duration (variation can be expected on time scale of the or- 
der of ro/c, i.e., C(lOms)), the black body spe ctrum is often 
sme ared. Even more than that, as was shown in lPe'er I (120081 ) 
and lPe'er fc Rvde l(|201H ). at any given instance, an observer 
sees simultaneously thermal photons emitted from different 
radii and different angles to the line of sight, hence having 
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different Doppler shifts. As a result, the expected photo- 
spheric emission is not a pure black body, but a combina- 
tion of black body spectra with different amplitudes. Thus, 
in fact, one expects to see a multi-color black body, as is in- 
deed seen. The full theory of multi-color black body emission 
from relativistically e xpanding plasmas recently appeared in 
iPe'er fc Rvdel l|201lh . 

The results of the fits to GRB090902B imply a very 
high value of e^, close to or even above equipartition, and 
an even higher value of e^, 85%- 95%. These values are 
much higher than the values predicted by the internal shocks 
model: the typical efSciency in energy dissipation by inter- 
nal shocks is no more than a few percent. The results obtain 
here (note that the exact nature of the dissipation process 
is not specified) thus raise another issue as to the valid- 
ity of the internal shock model scenario. This is added to 
GRB0 80916C, in which detailed analysis bv lZhang fc Pe'er I 
l|2009l ) concluded that an additional source of energy must 
exist between the photospheric radius and the dissipation ra- 
dius. The most plausible source of energy con sidered is mag- 
netic , i.e., a Poynting dominated outflow fsee lZhang fc Yan I 
I2OIII ). We can conclude that the high efficiency required, 
may hint toward Poynting dominated outflow in the case of 
GRB090902B as well, and may even be a general require- 
ment for all the bursts with pronounced high energy emis- 
sion observed by LAT. On the other hand, even if the outflow 
in GRB090902B is Poynting flux dominated, we do not ex- 
pect too high ratio of the Poynting to kinetic luminosity, a, 
since a ^ 1 results in suppression of the photospheric emis- 
sion, which is not observed in this burst. A more generalized 
treatment of photospheric models with arbitrary magnetiza- 
tion is outside the scope of this manuscript, and is left for 
future work. 

The separation made in this work into two emission 
zones, namely thermal emission originating from the pho- 
tosphere, and non-thermal emission originating from en- 
ergy dissipation at larger radii, provides a natural expla- 
nation to the delay o f the high energy photons, often seen 
in Fermi-LAT burs t s (lAbdo et al . ' 2009a: 'Rvde et aZ.1l2010l : 
iToma et all I2OI0I : iGhisellini et al. 2010 ). In our model, 
the non-thermal photons originate from dissipation above 
the photosphere, hence they are naturally seen at a delay 
with respect to the photospheric (thermal) photons, which 
are always the first to be observed. A pronounced thermal 
component at early times also provides a natural explana- 
tion to the harder slope at low energy observed during the 
first 1-2 s in many long bursts (jGhirlanda et al. Il2009h . In 
this general framework, the origin of the high energy, non- 
thermal photons can a l so be hadronic, as r ecently suggested 
(jRazzague et al. ll2010l : lAsano et al. II2OI0I ). Nonetheless, we 
showed here that a leptonic origin is consistent with the 
data. Our model has the advantage that it does not re- 
quire a large amount of energy in the hadronic component. 
Moreover, we did not specify the origin of the dissipation 
that lead to the emission of the high energy photons. Re- 
cently, there were sev eral suggestions of external shock ori- 
gin of these p hotons (Kum ar fc Barniol-Duran II2009I . |201C| : 
IGhisellini et a l. 2010). However, by fitting the broad band 
data we showed that excellent fits can be obtained if both 
the characteristic synchrotron breaks are below the Fermi 
energy band, i.e., less than ~ lOkeV. Thus, our model is 
consistent with internal dissipation origin of the non-thermal 



spectrum (see furt her discussion in IPiran fc Nakar I I2OI0I : 
IZhang et ar\\201'\\ ). 
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APPENDIX A: DECOMPOSITION OF THE 
SPECTRUM INTO ITS BASIC PHYSICAL 
INGREDIENTS 

The observed spectrum results from synchrotron emission, 
SSC and thermalization of the thermal photons, with a high 
energy cutoff resulting from pair production. As discussed in 
tj3.2l the three emission processes are expected to have simi- 
lar contributions to the observed spectrum of GRB090902B. 
Therefore, a "clean" decomp osition into the spec tral ingre- 
dients, as is presented in, e.g.. lSari fc Esin I (|200l|), does not 
exist in practice. 

The high energy spectral slope (above the thermal peak) 
does not have a direct correspondence to the power law in- 
dex of the accelerated electrons. Therefore, a decomposition 
of the spectrum into the basic radiative processes can be 
useful in understanding the physical processes that shape 
the spectrum. Fortunately, such a decomposition can (up to 
some level) be done numerically. 

The numerical results are presented in Figure lAll for 
dissipation radius r^ — 10^^ cm, and power law index p = 
2.2. The dash-dotted (red) curve represents the spectrum 
that would have been obtained if only synchrotron emission 
was considered. As the cooling frequency is low, the spectral 
index in the entire Fermi range is uF^, oc u^~^'^ oc r/""'^. 
When SSC is added (dashed, green curve) the spectral slope 
at low energies (below « 10* eV) is not affected, while the 
spectrum at high energies becomes nearly flat. This results 
from a combination of a decreasing synchrotron flux and an 
increasing SSC flux. Note that in this scenario, the flux at 
low energies is weaker than the flux obtained for the pure 
synchrotron scenario, due to the fact that part of the elec- 
tron energy is converted to SSC. Finally, the solid (blue) 
curve shows the combined effects when all the physical in- 
gredients are added. Inclusion of the thermal photons does 
not affect the low energy spectral slope (below the ther- 
mal peak). However, Comptonization of the thermal photons 
contribute to the more pronounced spectrum at high ener- 
gies. In addition, the inclusion of thermal photons leads to a 
sharper high energy cutoff, resulting from pair production. 
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